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MEASUREMENT  OF  THE  SHEAR  STRESS  ON  THE 
UNDERSIDE  OF  SIMULATED  ICE  COVERS 


Darryl  J.  Calkins  and  Andreas  Muller 


INTRODUCTION 

The  major  driving  force  exerted  on  an  ice  jam  is 
often  the  fluid  shear  transmitted  to  the  underside  of 
the  ice  cover.  Pariset  et  al.  (1966)  assumed  that  the 
ice  cover  shear  stress  was  equivalent  to  the  bed  shear 
stress  in  the  development  of  their  ice  cover  formation 
and  evolution  model. 

There  have  been  many  attempts  to  calculate  a  co¬ 
efficient  of  roughness  for  the  underside  o  fan  ice 
cover,  both  from  laboratory  and  field  observations 
(Nc/hikhovskiy  1964,  Ashton  ct  al.  1970,  Carey  1966 
and  1967,  Tatinclaux  and  Cheng  1978,  Ohashi  and 
Hanad  1970,  Tcsaker  1970,  Ismail  1972).  The  greatest 
difficulty  lies  in  measuring  the  physical  dimensions  of 
the  underside  roughness  element  of  the  ice  as  it  relates 
to  a  coefficient  of  friction. 

The  objective  of  the  research  described  in  this  re¬ 
port  was  to  experimentally  measure  the  drag  force 
and  hence  the  shear  stress  on  the  underside  of  simu¬ 
lated  ice  covers  in  a  laboratory  flume  using  square 
polyethylene  blocks  of  3,18-mm  thickness.  Velocity 
profiles  were  measured  beneath  the  fragmented  cover, 
and  the  shear  stress  was  calculated  from  tiicsc  profiles 
and  compared  against  the  measured  values. 


EXPERIMENTAL  APPARATUS 

The  experimental  work  was  carried  out  at  the 
Iowa  Institute  of  Hydraulic  Research  in  a  12.2-m-long, 
0.61  m-wide,  and  0.3-m-deep  tilting  flume,  maintained 
practically  horizontally  for  all  tests.  A  0.57-m-widc 
and  2. 44-m-long  wooden  frame  (Fig.  1 ),  constructed 
from  conventional  2X4  lumber,  was  suspended  3.05 
m  from  the  ceiling.  A  wire  screen  of  12.7-mm  mesh 
was  stapled  around  the  outer  surface  of  the  frame. 


excluding  the  upstream  end,  and  extended  approx¬ 
imately  0.1  m  below  its  base.  The  frame  was  designed 
so  that  the  mesh  would  not  come  in  contact  with  the 
flume  walls. 

Frame  displacement  was  measured  by  strain  gauge 
force  transducers  with  a  maximum  displacement  of 
0.06  mm.  The  two  transducers  were  connected  indi¬ 
vidually  by  thin  nylon  line  to  the  outside  members  of 
the  frame  at  the  upstream  section.  Initially  the  trans¬ 
ducers  were  calibrated  individually,  but  the  nylon  line 
stretched  so  that  the  system  as  a  whole  had  to  be  cali¬ 
brated.  Output  from  the  transducers  went  to  a  strip 
chart  recorder.  The  force  Ff  in  the  streamwisc  direc¬ 
tion  due  to  the  weight  of  the  frame  was  incorporated 
into  the  calibration  and  its  magnitude  can  be  neglected 

The  water  velocity  profiles  were  obtained  from  a 
Nova-Nixon  small  propeller  current  meter.  The  frag¬ 
mented  ice  cover  was  simulated  by  3.14-mm-thick 
X  31.7- X  38.1 -mm  polyethylene  high  density  plastic 
blocks  of  0.92  specific  gravity. 

EXPERIMENTAL  PROCEDURES 

A  constant  water  level  of  0.152  m  was  maintained 
tor  all  experiments  and  after  some  preliminary  tests 
two  flow  velocities  were  predominantly  used,  0.082 
and  0.1  3  m/s.  The  critical  velocity  for  submergence 
of  an  individual  block  was  calculated  to  be  approxi¬ 
mately  0.13  m/s  based  on  data  from  Larsen  (1975) 
and  Tatinclaux  ct  al.  (1977). 

When  the  desired  flow  velocity  was  obtained  in  the 
f  lumc,  the  nylon  lines  between  the  transducers  and 
frame  were  brought  into  tension  by  moving  the  pair 
of  transducers  in  an  upstream  direction  until  a  small 
force  registered  on  the  recorder  was  zeroed  out.  Ice 
was  introduced  into  the  channel  near  the  headbox  from 
a  chute  angled  from  the  support  rails  to  the  water  sur- 


Flow 


Figure  1.  Fxperimentul  setup  in  flume. 


surface.  When  the  ice  dropped  into  the  fow,  surface 
waves  were  propagated  downstream,  setting  up  temp¬ 
orary  oscillations  in  the  frame  as  indicated  by  the 
recorder  output. 

The  total  force  acting  on  the  ice  mass  used  for  cal¬ 
culations  was  determined  when  the  system  had  stabil¬ 
ized  and  a  constant  reading  attained.  Velocity  profiles 
were  measured  0.76  m  upstream  from  the  downstream 
edge  of  the  ice  accumulation.  The  weight  of  ice  dis¬ 
charged  into  the  flume,  length  of  accumulation  and 
averagejhickness  were  also  measured.  The  mean  ve- 
loeity  Vm  in  the  channel  was  calculated  from  the  dis¬ 
charge  measurements  obtained  with  an  in-line  orifice 
plate. 

ANALYSIS  OF  FORCES 

The  ice-collecting  frame  was  designed  to  exclude  the 
friction  force  that  could  be  developed  on  the  walls  of 
the  flume  by  ice  pieces.  The  total  force  FM  measured 
by  the  transducers  in  the  direction  of  flow  consisted 
of  three  components:  1 )  hydrodynamic  or  form  drag 
Fh  on  the  leading  edge,  2)  shear  force  Fs  due  to  the 
moving  fluid  beneath  the  cover  including  form  drag 
of  the  ice  blocks,  and  3)  weight  of  the  ice  mass  f  ], 
f  igure  I  shows  a  schematic  view  of  the  frame  and  the 
forces  acting  on  the  system. 

The  forces  are: 

^  1  F\  ( 1 ) 

where  F"H  =  Vip,  Cd  t  B  (hydrodynamic  force) 

Fj  =  ryLB  (shear  force) 

F | |  =  H’|  sin (3  (cover  force) 

P\  -  density  of  fluid 
Cd  =  drag  coefficient  =  1 .0 
r  ~  thickness  of  ice 
B  =  width  of  frame 
L  -  length  of  ice  accumulation 


=  mean  fluid  velocity 

r,  =  shear  stress  on  the  ice  cover  underside 

The  slope  of  the  water  surface  gradient  was  approx¬ 
imately  5v  If)-'  m/m,  and  the  force  F{  in  the  down¬ 
stream  direction  due  to  an  ice  cover  was  small  1  mN) 
and  consequently  could  be  excluded.  It  will  be  shown 
in  the  experimental  data  section  that  the  hydrodynamic 
force  Fn  is  small  compared  to  the  measured  total  forces 
except  possibly  lor  short  cover  lengths.  Therefore,  the 
total  force  recorded  In  the  transducer  consisted  pri¬ 
marily  ol  the  shear  stress  transferred  to  the  cover's 
underside. 

The  calculation  ot  the  shear  stress  on  the  ice  cm  -  ’s 
underside  was  based  on  the  projected  planar  area  of  the 
cover  within  the  ice  frame,  such  that 

l/fi  (2) 

where  the  subscript  m  refers  to  the  measured  value. 

The  measured  shear  stiess  was  compared  with  sheer 
stress  calculated  Irom  the  velocity  profile  on  the  assump¬ 
tion  that  the  logarithmic  velocity  distribution  could  be 
applied  to  a  How  depth  only  a  few  times  the  si/e  of 
the  roughness  scale.  A  further  assumption  was  made 
that  the  flow  cross  section  under  the  ice  cover  could 
be  divided  into  two  /ones  Irom  the  plane  ot  maximum 
velocity  I )  the  ice  /one  and  2)  the  bed  /one.  Based 
on  this  division  of  areas,  the  logarithmic  velocity  dis¬ 
tribution  was  applied  to  each  section  individually, 
f  he  only  difficulty  lay  in  assessing  the  position  of 
the  upper  boundary  since  the  fragmented  cover  had 
such  a  variable  thickness. 

EXPERIMENTAL  RESULTS 

An  initial  set  ot  experiments  was  conducted  to 
determine  the  necessary  cover  length  that  would  en¬ 
sure  a  proper  boundary  layer  development  and  to  in- 
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Table  1.  Initial  shear  stress  tests. 


Run 

l  ertifth  ot 
ice  on  er  l. 

( ni  1 

H't'ujht  ot 
cow  W 
tM 

/  tiu  kne s.% 

of  core/  l 

(mm ) 

rota! 
tone  1  M 
(mN) 

Mean 

edacity  Vm 
(m/s) 

Measured 

shear  stress  r  - 
(Mm7) 

Normalized 
thickness  t/v  n 

friction 
factor  fjm 

Note 

SI 

O.o  1 

8.89 

3.18 

34.2 

0.1 19 

0.0948 

. 

S2 

1.22 

13.34 

3.18 

46.7 

0.1 19 

0.0868 

-• 

- 

S3 

1.83 

20.09 

3.18 

8  3.6 

0.119 

0.077  3 

- 

- 

Si 

2.44 

35.59 

3.18 

106.3 

0.119 

0.0737 

- 

- 

S' 

1.58 

35.59 

(13.3) 

99.2 

0.082 

0.1058 

0.095 

0.1259 

1 

So 

1 .92 

48.93 

05.1) 

124.1 

0.082 

0.1092 

0.110 

0.1299 

i 

S7 

2.*S 

62.28 

(15.5) 

152.1 

0.082 

0.1082 

0.1 14 

0.1287 

1 

ss 

2.44 

35.59 

3.18 

233.5 

0.137 

0.1618 

- 

- 

S9 

2.44 

35.59 

3.18 

193,5 

0.137 

0. 1  34  1 

- 

4 

S 10 

2.44 

35.59 

3.  IS 

164.1 

0.137 

0.1  140 

- 

5 

SI  1 

2.44 

*5.59 

3.18 

134.8 

0.137 

0.0934 

- 

- 

si : 

2.44 

35.59 

3.18 

134.8 

0.137 

0.0934 

- 

- 

6 

SI  t 

2.44 

35.59 

3.18 

1  14.8 

0.137 

0.0934 

- 

- 

7 

S14 

2.1  1 

5  3.38 

(14.7) 

289.1 

O.i  19 

0.2296 

0.106 

0.1297 

1 

SIS 

2.15 

5  3.38 

(14.7) 

321.6 

0.119 

0.2555 

0.106 

0.1443 

1 

S  1 6 

2.2*> 

5  3.38 

(1  3.0) 

105.4 

0.079 

0.0778 

0.098 

0.1010 

1 

si: 

2.29 

1.2.28 

(15.9) 

148.3 

0.079 

0. 1 096 

0.114 

0.1405 

1 

sis 

2.29 

(.2.28 

(25.9) 

140.1 

0.079 

0.1035 

0.114 

0.1327 

1 

SI 

2.41 

106.7b 

(  8.5) 

676.1 

0.119 

0.4747 

0.204 

0.2681 

1 

S20 

2.44 

35.59 

l  8.5) 

163.2 

0.137 

0.1132 

0.059 

0.0482 

S2I 

2.44 

35.59 

3.  IS 

11  9.2 

0.137 

0.0827 

0.059 

0.0352 

8 

LI 

0.0  J 

10.32 

3.18 

1  1.5 

0.079 

0.0319 

- 

- 

l  2 

1.22 

20.04 

3.18 

14.2 

0.076 

0.0197 

- 

It 

1.81 

30.0b 

1.18 

22.7 

0.079 

0.0210 

- 

~ 

L4 

2.44 

■11.28 

3.18 

27.0 

0.079 

0.0187 

- 

1.5 

2.44  + 

.1.18 

28.3 

0.079 

0.0196 

- 

2 

L6 

2.44  + 

3.18 

28.3 

0.079 

0.0196 

- 

3 

t.  7 

2.44 

3.J8 

83.2 

0.1  19 

0.05  7  7 

- 

3 

L8 

2.44 

3.18 

83.2 

0.1  19 

0.0577 

- 

- 

3 

L9 

2.44 

3.18 

J  54.8 

0.140 

0.1074 

3 

1  II) 

2.44 

41.28 

1.18 

176.6 

0.168 

0.1225 

-• 

3 

LI  1 

0.0 1 

10.32 

3.18 

12.9 

0.082 

0.0358 

- 

- 

112 

1.22 

20.  h4 

1.18 

1 5.6 

0.082 

0.0216 

1. 1  t 

1.83 

30.9b 

3.18 

22.7 

0.082 

0.0210 

- 

Notes:  !.  I  hickness  iorrevted  bv  porosils  -  0.67. 

2.  (‘over  length  inuejsed  hv  0.0|  in  upstream  ot  frame  (0.91  xO.6  I  -m  sheets). 

3.  Cover  length  increased  hv  1.82  nt  upstream  of  (tame  (0.91  xO.  6  l-m  sheets!. 
■1.  t  ive  blocks  removed  from  downstream  screen  lace  in  run  S8. 

s.  I  en  blocks  removed  Corn  downstream  screen  lace  in  run  98. 
b.  (  over  length  increased  b\  0.91  m  upstream  (single  blocks). 

7.  Cover  length  increased  hv  1.82  rn  upstream  (single  blocks). 

8.  three  blocks  at  45  to  60’  removed  trom  run  S20, 


vestigate  it  one  can  neglect  the  hydrodynamic  force 
on  the  leading  edge.  I  he  data  for  these  experiments 
can  be  found  in  Tables  I  and  2.  Both  large  ice  sheets 
(O.bl-x  0.59-m-widc)  and  ice  blocks  (all  3.18  mm  in 
thickness)  were  tested. 

From  the  initial  experiment  it  was  determined  that 
a  cover  length  greater  than  2  m  was  necessary  to  have 
a  proper  boundary  development  and  that  both  the 
hydrodynamic  force  at  the  leading  edge  and  the  weight 
of  the  ice  in  the  downstream  direction  can  he  excluded 


because  of  their  minor  contributions  compared  to  the 
shear  force. 

As  a  check  on  the  sensitivity  of  the  experimental 
apparatus,  smooth  flat  sheets  of  0.1 6-x  0.59-m-wide 
polyethylene  plastic  were  placed  in  the  free-swinging 
frame  and  the  drag  force  was  measured  for  two  flow 
velocities.  In  a  discussion  of  the  turbulent  boundary 
layer  over  smooth  flat  plates,  Schlichting  (1968)  de¬ 
rived  the  skin  friction  coefficient  based  upon  the  log¬ 
arithmic  velocity  distribution  and  the  drag  force  acting 
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Table  3.  Measured  and  computed  drag  forces  on  the 
smooth  ice  sheets. 
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on  a  plate  ot  /eio  incidence,  However,  the  Reynolds 
numhet  In;  the  plate  lengths  and  velocities  considered 
lies  in  the  tiansiiimi  /one.  and  one  must  determine 
whether  the  hound.nv  laser  is  laminar,  turbolent  ot  in 
transition  over  the  plate  length. 

It  was  iniiiailv  assumed  that  the  llovc  was  probable 
turbulent  because  the  icc  sheets  had  blunt  square  lead¬ 
ing  edges,  therein  tupping  the  houndat \  laser  into  the 
turbulent  ,,eginie,  and  this  assumption  was  later 
con  I  n  med. 

Sell  lie  111  trig  i  I  9l>8)  I  it  ted  l  he  relation  between  C. 

I  tie  toelfic  rent  t»!  skin  trie  I  ion  and  :'t  the  Res  Holds 
number  as  developed  horn  the  logarithmic'  uTocits 
protrle  into  a  more  convenient  empn  real  lorril  tor 
turbulent  l lows  suJ:  that 

c,  0.455  (log .  ..  A*  )- os  (3) 

where  1 1  /  I  !■  arid  r  is  the  kinematic  vrscosits . 

lo  compare  the  measured  and  calculated  values  ot 
the  total  drag  acting  on  a  smooth  plate,  tests  V  I  and 
V2  we>e  pet  lonnc'd  (see  I  tbie 2)  I  he  drag  on  the 
snioc.h  sheets  / i  is  given  hs 

p,  t  ■  I  HI.  (4) 

I  able  I  gives  the  tneasnifd  drag  lorceand  the  com¬ 
puted  values  1 1 >r  these  two  test  mils. 

t  he  measured  values  ol  the  drag  lorae  are  in  tela- 
lisels  close  agreement  with  those  values  obtained  using 
the  Prandti-Schlic  ting  turbulent  flow  relationship  lor 
the  skirr  Ir  ictr  >rr.  Based  on  this  into/ illation,  it  was 
conrluded  iliat  the  apparatus  was  capable  ol  measur¬ 
ing  quite  accuratelv  the  shear  stress  exerted  In  the 
lluid  to  a  iloating  cover  and  that  the  boundary  laser 
was  in  lad  t nr  Indent. 

1  he  el  tec  l  .4  ice  cos  cl  length  in  the  trame  on  the 
drag  lorce  measured  is  indirectly  shown  in  I  igure  2. 

I  or  the  1 1  -w  selocils  lests  with  short  cover  lengths 
the  measured  torces  arc  extremely  low  (  15  mNj.and 

the  et tec  t  oi  the  IrvdrocK  narnic  lor  cc  and  weight  ot 
the  ice  cover  Incomes  nicieasinglv  important  with 
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/  inure  2.  I  he  decreasing  effect  ot  the  lead- 

inn  edge  form  drag  with  increasing  cover 

length  as  measured  by  the  overall  shear  stress. 

decreasing  length.  I  or  example,  1. 1  with  a  cover 
length  of  0.6 1  niN  has  a  potential  conn  ihution  from 
the  hydrodynamic  force  of  roughly  6  mN,  using  c()  = 

0.1 .  Consequently,  short  cover  lengths  would  give 
higher  values  ol  shear  stress  il  the  hydrodynamic  force 
were  not  subtracted,  and  tins  is  shown  in  I  igure  2.  The 
value  ot  the  effective  shear  sir  ess  decreases  with  increas¬ 
ing  cover  length  for  tests  1. 1  and  12  but  remains  nearly 
constant  lor  L3-I.6. 

I  or  runs  L5  and  Lb,  additional  sheets  were  placed 
ininiediatcTs  upstream  ol  the  Iramc,  such  that  a  proper 
boundary  laser  and  velocity  distribution  might  be 
achieved.  I  hese  additional  sheets  did  not  come  in  con¬ 
tact  w  ith  the  suspended  Iramc,  and  lor  the  low  veloc- 
its  I  low  under  the  large  sheets  the  shear  stress  was 
neatly  constant  aftet  the  cover  had  reached  a  length 
ol  rouglds  2  m.  In  experiments  SI  1-S13  (small  blocks 
ol  one  laser  thickness),  the  covei  was  extended  up¬ 
stream  from  the  trame  approximately  0.9  and  1.8  m 
respectively ,  and  again  a  constant  value  of  shear  stress 
was  maintained.  The  data  lor  the  higher  velocity  flow  s 
(SI -54  and  S II -SI  3)  seem  to  suggest  that  a  cover 
length  ol  approximately  2.4  m  was  satisfactory  lor 
measuring  the  shear  stress  and  the  ettects  ot  the  hydro- 
c> narnic  force  are  now  small  compared  to  the  shear 
lorce. 

I  he  velocity  proliles  at  0.76  m  I  mm  the  downstream 
end  o.  the  Iramc  are  shown  in  I  igure  3  for  increasing 
cover  lengths  upstream,  these  profiles  reveal  that,  at 

higher  velocities  (Ml. 1 3  m/s),  a  cover  length  of  2.5  to 
3.5  m  is  necessais  for  proper  boundary  layer  develop¬ 
ment  at  the  velocity  measuring  position. 

A  relative  minimum  thickness  ot  the  ice  jam  can  be 
computed  on  the  basis  ol  the  surlace  area  ol  icc  in  the 
frame  and  the  total  weight  introduced  upstream.  This 
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Figure  3.  Velocity  profiles  tit  0.  76  m  from  the  downstream  end  for 
increasing  cover  lengths  upstream. 


would  represent  the  minimum  thickness  obtainable 
with  no  allowance  for  porosity.  It  will  be  shown 
that  the  porosity  remains  relatively  constant  tor  all 
tests.  The  absolute  roughness  on  the  underside  of  the 
accumulated  ice  was  not  measured  because  time  did 
not  permit  extensive  measurements  during  these  tests. 

The  orientation  of  the  blocks  in  the  ice  cover  at 
the  downstream  end  of  the  frame  had  a  significant 
effect  on  the  total  force  measured.  Tor  example,  run 
S20  had  three  blocks  against  the  downstream  wire 
screen  at  45°-60°;  removing  only  these  three  blocks 
resulted  in  a  27%  decrease  in  the  force  exerted  on  the 
system  as  measured  in  run  S21 .  This  same  observation 
was  noted  in  run  S8,  where  five  blocks  were  removed 
at  the  wire  screen  for  run  S9  and  five  more  for  run 
S10,  effecting  a  30%  reduction  in  total  force  from  run 
S8. 

For  run  S20,  the  computed  drag  force 

Ft  --;(pV*ACd)l2  (5) 

on  these  three  blocks  at  60°  is  approximately  30  mN 


(where  Cd  =  1 .0  and  A  is  the  projected  area  of  the 
block).  Comparison  of  this  value  with  the  loss  in  force 
of  44  mN  measured  by  the  system  when  the  blocks 
were  removed  indicates  that  a  few  blocks  can  account 
lor  a  major  portion  of  the  force  due  to  large  areal 
projections.  These  two  observations  clearly  show  that 
care  must  be  taken  at  the  downstream  end  so  that 
unobstructed  ice  floes  do  not  project  excessively  into 
the  flow  at  steep  angles  and  create  excessive  drag 
forces  compared  to  the  skim  drag  of  the  accumulation 
thickness. 

For  short  ice  cover  lengths  3.0),  the  hydro- 
dynamic  drag  Fh  appears  to  be  an  important  contrib¬ 
utor,  but  as  the  cover  length  increases,  the  ratio  of  the 
hydrodynamic  form  drag  on  the  leading  edge  to  the 
total  load  is  small  (<  0.05).  This  is  because  the  ice 
cover  underside  drag  force  is  significantly  greater  than 
the  ice  cover  frontal  form  drag. 

The  velocity  profiles  beneath  the  simulated  ice 
covers  are  quite  similar  except  that  the  roughness  of 
the  respective  boundaries  influences  the  position  of 
the  maximum  velocity.  Figure  4  shows  the  velocity 
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/  iijurr  i  \  <7i>(  it\  pinhlcs  at  0.  '(>  m  Iron)  downstream  end  be¬ 
neath  :<uh/h  n  (• ,  in  i' s. 


profiles  .1  ml  the  lonesponding  change  m  position  o| 
t he  maximum  velot il\  due  to  hull)  inci  casing  nun 
thickness  .iikI  boundarv  roughness.  Because  ihr  Ilium- 
bottom  was  relatively  smooth,  the  position  ol  ma\i 
mum  velocilv  with  llit-  latgc  ice  sheets  was  neaih  one 
hall  the  depth  ol  Mow  beneath  thecovet.  As  the  cover 
mereased  in  toughness  (eteated  h\  mtrodiH  mg  greater 
quant  ities  ol  ne  blocks  pet  unit  time  upstream),  the 
position  ol  ma Mimim  velocits  approached  the  smoother 
boundary 

ANALYSIS  OF  DATA 

I  he  data  m  I  allies  I  and  2  teptesent  the  results  ol 
i!ie  shear  stu-ss  measurements  on  single  thickness 
large  sheets  (I  ),  small  blocks  (S)  and  on  both  the  sheets 
and  blocks  when  the  velocits  protile  measurements 
were  taken  (V). 

I  ,  is  the  mean  velocits  ol  the  Now  beneath  the  ice 
cover  in  the  region  liom  the  point  ol  masimum  veloc 


its  to  the  bottom  ol  the  ice  cover  as  determined  Irom 
the  velocits  pioliles.  I  lie  shear  velocits  in  the  ice  sec¬ 
tion  ’  was  determined  from  the  velocity  profiles,  as¬ 
suming  a  logarithmic  velocits  distribution  Irons  the  ice 
bottom  to  the  point  ot  maximum  velocity. 

I  he  biggest  difficulty  in  calculating  the  shear  stress 
lies  in  assessing  the  position  of  zero  velocity  with  re¬ 
spect  to  the  boundary  of  the  underside  of  the  ice 
cover.  Since  the  boundary  is  very  irregular,  it  is  ex¬ 
tremely  difficult  to  locate  from  measurements.  Fur¬ 
thermore,  a  question  arises  as  to  where  one  should 
take  the  velocity  profile  measurements.  Hence,  the 
calculation  ol  the  friction  velocity  from  the  velocity 
profiles  is  critical,  as  the  position  ot  the  boundary  can 
significantly  influence  the  slope  of  the  velocity  pro¬ 
file  gradient.  The  shear  stress  computed  from  the  re¬ 
lation  rj(,  u , * uses  the  friction  velocity  Uj*  as 
calculated  Irom  the  velocity  prolile  gradient. 

Ihe  general  equation  for  Darcy's  coefficient  of 
(r  id  ion  t  is 
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Figure  5.  Definition  sketch  of  flow  symbols. 


Figure  6.  Relationship  between  computed  and  measured  values  of  the 
ice  cover  shear  stress. 


f-  8tIp(V^.  (6) 

The  specific  resistances  f[t  are  calculated  using  the  mean 
velocity  in  the  ice  section  V  v  and  fm  is  the  calculated 
friction  coefficient  using  the  mean  channel  velocity 
Vm  (see  Fig.  5  for  illustration  of  terms). 

The  porosity  p  of  the  ice  accumulation  was  calcu¬ 
lated  from  the  data  contained  in  Table  2.  The  icc 
cover  porosity  is  defined  as  the  ratio  of  the  volume  of 
the  voids  filled  with  water  to  the  total  volume.  The 
porosity  values  arc  consistent  with  those  obtained  by 
Tatinclaux  and  Cheng  (1978). 

Figure  6  shows  the  comparison  of  the  shear  stress 
calculated  from  eq  2  and  that  computed  from  the  ve¬ 
locity  profile  reduction.  For  the  low  velocity  runs 
(0.082  m/s)  the  agreement  between  the  measured  and 
calculated  shear  stress  was  adequate,  but  the  higher 
velocity  tests  showed  an  increasing  divergence  of  al¬ 
most  2:1  for  computed  vs  measured  values. 


Several  investigations  have  shown  that  the  logarith¬ 
mic  distribution  fails  to  conform  to  the  observed  ve¬ 
locity  profile  when  large  discrete  roughness  elements 
a  e  present  at  the  boundary  (Davar  and  Ismail  1977, 
O'Laughlin  1965,  and  Ismail  1977).  Although  the  log¬ 
arithmic  velocity  equation  appeared  to  be  a  reasonable 
fit  of  the  velocity  profile,  the  von  Karman  constant 
k  (i.c.  the  ratio  of  friction  velocity  and  the  logarithmic 
gradient  of  the  mean  velocity)  seems  to  vary. 

There  arc  four  possible  explanations  for  the  dis¬ 
crepancy  in  the  shear  stress  data  at  the  higher  veloc¬ 
ity  tests: 

1.  The  propeller  current  is  not  small  compared  to 
the  flow  depth. 

2.  The  propeller  many  not  be  responding  properly 
to  shear  flow. 

3.  The  porous  boundary  may  be  important  and 
three  dimensional  flow  may  be  present. 

4.  The  exact  position  ot  the  zero  velocity  near  the 
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/  /(jure  7.  typical  tril  lion  factor  relationship 
w  ith  normalized  ice  cover  thickness  using  V,„. 


f  igure  8.  typical  friction  coefficient  relationship 
with  normalized  ice  cover  thickness  using  V(. 


porous  boundary  is  in  question;  Chu  and  Gelhar 
(1972)  obtained  a  von  karman  constant  of  0.27 
above  a  porous  wall. 

Sufficient  information  was  obtained  to  calculate  the 
coefficient  ot  friction  factor  f  for  the  tests  based  on  the 
frame  force  measurements.  Two  values  of  the  friction 
factor  were  calculated  for  runs  VI  -V 1 2,  one  using  the 
mean  upstream  open  water  flow  velocity  f  m  and  the  sec¬ 
ond  the  mean  velocity  ol  the  flow  in  the  ice  region  V.. 

Normalizing  the  measured  thickness  by  the  upstream 
depth,  two  plots  of  /|m  and  /• f  vs  t/y ■  are  presented  (Tig. 

7  and  8),  which  generally  confirm  the  assumption  that 
the  roughness  ol  the  cover  increases  with  increasing 
cover  thickness.  The  mean  velocity  in  the  ice  section 
was  usually  lower  than  the  upstream  open  flow  mean 
velocity  which  leads  to  slightly  higher  friction  factors. 

I.arlier  data  were  collected  on  the  ice  cover  shear 
stress  where  the  thickness  of  the  cover  was  not  meas¬ 
ured,  and  an  attempt  was  made  to  recover  the  estimated 
thickness.  In  runs  V5-V 1 2,  excluding  V9,  the  porosity 
of  the  accumulation  an  be  calculated  from  measure¬ 
ments  ol  the  average  ice  thickness  and  volume  of  ice 
in  the  test.  Taking  the  average  porosity  value  of 
p  0.67,  the  expected  thickness  is  shown  in  parenthe¬ 
ses  in  Table  I . 

On  the  assumption  that  the  roughness  increases  in 
proportion  to  the  ice  cover  thickness  for  these  data, 
the  friction  factor  f-  can  be  plotted  vs  the  nondimen- 
sional  roughness  parameter  t!yn  as  is  shown  in  figure 
9.  Using  all  the  data  from  Tables  1  and  2  with  accum¬ 
ulation  lengths  greater  than  2  m  and  with  more  than  a 
single  thickness  ice  cover,  it  can  be  seen  that  the  data 
lie  within  the  same  general  zone  as  that  shown  by 


Tatinclaux  and  Cheng  (1978).  A  regression  lit  on  the 
new  data  was  omitted  because  they  appear  to  fall  with¬ 
in  the  same  general  range  as  those  of  previous  investi¬ 
gators,  although  maybe  slightly  lower. 
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figure  9.  Relationship  between  Darcy’s  f  and 
the  normalized  ice  cover  thickness  compared 
with  latinclaux  and  Cheng  's  1978  data  points 
which  are  represented  by  a  shaded  area  con¬ 
taining  90%  of  their  data. 
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/  igure  1 0.  The  friction  far  tor- Reynolds  number  re¬ 
lationship  for  open  channels  with  different  flow  re¬ 
gimes  (see  Chow  1959). 

I  he  assumption  that  the  roughness  of  the  under¬ 
side  of  the  fragmented  ice  cover  increases  with  the 
cover  thickness  is  not  a  strong  indicator  of  the  rough¬ 
ness  but  appears  adequate  lor  this  single  block  si/e  ex¬ 
periment.  f  or  all  future  experiments,  the  roughness 
heights  should  be  measured  and  used  in  place  of  the 
ice  cover  thickness,  and  field  measurements  of  ice  jam 
roughness  elements  would  be  extremely  valuable  to 
guide  furthei  work. 

The  friction  (actor  for  the  various  types  of  ice  cover 
I )  large  sheets  2)  single  thickness  blocks  and  3)  mul¬ 
tiple  thickness  covers  are  plotted  vs  the  Reynolds  num 
her  Rt  -  VJ/v,  where  /  is  a  characteristic  length  chosen 
to  be  the  hydraulic  radius  R  in  Figure  10  (from  Chow 
1959). 


The  data  lor  the  large  smooth  sheets  lie  near  the 
Blasius  equation  for  turbulent  flow  over  flat  plates, 
while  a  cover  ol  a  single  thickness  of  small  blocks  shows 
a  slightly  higher  resistance.  The  scatter  in  the  multiple 
floe  thickness  data  lies  in  the  regions  of  the  various 
Ba/in  channels. 

Ne/hikhovskiy  (1964)  presented  data  on  the  Man¬ 
nings  roughness  coefficient  n  and  thickness  of  ice 
sheets  for  rivers.  His  data  points  are  replotted  in  Fig¬ 
ure  1 1  along  with  the  laboratory  data  of  this  study. 
Both  sets  ol  data  follow  the  general  trend  of  increasing 
roughness  coefficients  for  greater  ice  block  thicknesses. 
It  is  clear  that  the  two  data  sets  are  in  different  regimes, 
but  the  roughness  coefficients  increase  at  nearly  the 
same  exponential  tate  of  0.4. 

CONCLUSIONS 

The  method  presented  for  measuring  the  shear  stress 
on  the  underside  of  fragmented  covers  in  small-width 
channels  eliminates  the  effect  of  wall  friction  force 
which  is  significant.  The  shear  stress  computed  from 
the  velocity  profiles  for  the  higher  velocity  experiments 
was  on  the  order  of  twice  the  measured  values,  indi¬ 
cating  a  deficiency  in  the  von  Karman-Prandtl  logarith¬ 
mic  velocity  distribution.  This  effect  may  stem  from 
1)  inability  to  adjust  the  profile  to  the  boundary  where 
the  velocity  is  zero,  2)  three-dimensional  flows  around 
the  large  roughness  elements,  or  3)  the  possibility  of 
the  porous  wall  influence. 

The  nondimcnsional  friction  factor  /  based  on  the 
thickness  of  the  fragmented  ice  cover  was  adequate  in 
assessing  the  roughness,  indicating  that  the  roughness 
was  proportional  to  the  increasing  thickness  of  the 
cover.  The  values  for  the  friction  coefficient  of  the 
model  ice  arc  in  agreement  with  various  reported 
values  for  real  ice  covers  (Carey  1966,  1967,  Tcsaker 
1970,  Nezhikhovskiy  1964). 


Ice  Thickness  (m) 


Figure  / 1.  Lmpirical  relationship  between  roughness  coefficient  and  thickness  of  the  ice  cover. 
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